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ABSTRACT 

Be stars are rapidly spinning B stars surrounded by an outflowing disc of gas in 
Keplerian rotation. Be star/X-ray binary systems contain a Be star and a neutron 
star. They are found to have non-zero eccentricities and there is evidence that some 
systems have a misalignment between the spin axis of the star and the spin axis 
of the binary orbit. The eccentricities in these systems are be caused by a kick to 
the neutron star during the supernova that formed it. Such kicks would also give 
rise to misalignments. In this paper we investigate the extent to which the same 
kick distribution can give rise to both the observed eccentricity distribution and the 
observed misalignments. We find that a Maxwellian distribution of velocity kicks with 
a low velocity dispersion, cr^ « 15kms~^, is consistent with the observed eccentricity 
distribution but is hard to reconcile with the observed misalignments, typically i ^ 
25°. Alternatively a higher velocity kick distribution, ak = 265kms~^, is consistent 
with the observed misalignments but not with the observed eccentricities, unless post- 
supernova circularisation of the binary orbits has taken place. We discuss briefly how 
this might be achieved. 

Key words: stars: emission lines, Be; stars: neutron; X-rays: binaries; accretion, 
accretion discs; supernovae 



1 INTRODUCTION 

. Be stars were discovered by ISecchil |l863) who observed 
' emission lines in 7 Gas. These stars are ra pidly rotatin g at 
about 70 per cent of their break up velocity (|Porteij| 19961 1. In 
fact, relative to their break up velocities, they are the fastest 
rotating bodies observed. They are early type main-sequence 
stars which have shown Ha in emission at least once. They 
are variable in brightness and spectra which show broad Hel 
absorption and emission at either visual or UV wavelengths. 
It is found that the emission, and so presumably the discs, 
are only temporary and so Be stars become B stars and vice 
versa. 

If the disc is viewed edge on the Be star is seen as a shell 
star. The spectra then show Balmer emission with sharp ab- 
sorption cores, narrow absorption lines of ionized metals and 
broad Hel absorption. Be-star dis cs vanish and re-appear on 
timesc ales of a few hundred days. lOachs. Kiehling fc Engela 
l|l98^ studied Balmer emission line profiles and concluded 
that the envelopes surrounding Be stars are in Keplerian mo- 
tion within the disc. The disc shows optical and IR emission 
lines and an IR continuum excess. 

The stars 7 Gas and 59 Gyg have shown two successive 
shell events. These were associated with a remarkably syn- 
chronous quasi-cyclic variation of the emission line width in 
all observed emission lines that has been called spectacular 



variation iHummelll 19981) . The change in emission line width 
removes the correlation between the projected surface veloc- 
ity, usini, and the FWHM (full width half maximum) and 
so a circumstellar equatorial disc fails to explain the spectac- 
ular variations. The emission lines and shell Be stars are ex- 
plained by differences in disc inclination to the lin e of sight, 
so transitions between the two were not expected. iHummell 
(1998) explains the spectacular variations by a Keplerian 
disc which is somehow tilted with respect to the equatorial 
plane of the star. The variation in emission line widths and 
profile shapes are then due to the precession of the disc. He 
suggests that the sequence of alternating shell-phases and 
narrow single peak phases is due to the variation in disc in- 
clination caused by precession. The idea that a disc might 
change its inclination to the line of sight is bo rne out by 
observations of 28 Tan (Pleione) bv lHiratal (|2007l ). This star 
also changes between B star. Be star and shell star and in 
this star the intrinsic polarization angle changes in phase 
with the se variations. Th e cause of the precession is also not 
clear but lHumme il (|1998|) suggested that it might be induced 
by tides from a binary companion. The two systems 7 Gas 
and 59 Gyg are binary. The system 28 Tau appears to show 
radial velocit y variations, although there is no confirmed or- 
bital period i|Rivinius. Stefl fc Baad3l2006l '). In both 7 Gas 
and 28 Tau the misalignment angle between the stellar equa- 
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tor and the disc/orbital plan e is thought to be around 25° 
l|Hummellll998l : lHiratall2007h . 

In further support of this possibiUty, we note that the B- 
star binary PSR J0045-7319 has a spin- orbit misalignment 
suRgested by its orbital plan e precession (|Kaspi et al.lll99(j : 
iLai. Bildsten fc Ka psi"l995'). This misalignment in a B-star 
orbit means that misalignment in Be stars in not uncommon. 
In this ca se the B st ar rotates retrogradely with respect to 
the orbit (|Lailll996al l. 

The standard model for Be-star discs is that they are 
decretion discs with the mass expelled from th e neighbour- 
hood of the Be star itself (|Cassinelli et al.ll2002h . In this case 
we expect the plane of the inner disc to be aligned with the 
spin axis of the Be star. There are two reasons for the disc to 
be fou nd at an angle discussed in the literature. First IPorterl 
l|l998h suggested that the disc warping and pr ecession migh t 
be caused by a radiation-induced instability (|Pringlell 1999 1. 
Secondly, as we reported above, it is widely suggested that 
the warp and precession are caused by a misalignment be- 
tween the spin axis of the Be star and the orbit of the binary 
companion. If the disc is a decretion disc we expect the inner 
edge to be aligned with the equatorial plane of the B star 
and the outer edge to tend to be tidally aligned with the 
orbital plane. Thus there must be a warp at some radius in 
the disc. Here we focus on this second possibility and con- 
centrate on the Be/X-ray binaries in which the companion 
stars are neutron stars. 

Neutron stars, observed as radio pulsars, have 
space velocities m uch greater than their progenitors 
l|Gunn fc Ostrikeil Il970l ). The accepted explanation 
for this is that supernova explosions are asymmetric 
and give ve ry large k i cks t o the new l y for med neu- 
tron stars (|Shklovskiil Il970l : ISutantvd Il978l ). Indeed 
some supernova rem nants show evidence for 
metric explosions |Morse. Winkler fc Kirshneij 
lAschenbach. Egger fc Trumpet 19951 ). Thus any 
that contains a neutron star could have had a supernova 
kick. It is also found that Be/X-ray binaries have high 
eccent ricities that cannot be explaine d without supernova 
kicks ijVerbunt fc van der Heuveilll995l l. 

In a Be-star binary system, prior to the supernova in 
which the core of its companion collapses to a neutron star, 
we expect the Be-star spin to aligned with the orbit and 
for the orbit to be circular. The kick from the supernova 
has two effects, it makes the orbit eccentric, and perhaps 
even unbinds it, and it misaligns the orbit with the spin 
axis of the Be star. Thus information about the distribution 
of eccentricities in Be-star systems can in principle give us 
information about the distribution of spin-orbit misalign- 
ments. 

iLvne fc Lorimeij (| 19941 1 analysed the known pulsar ve- 
locities and concluded tha t they were born w i th a mean 
speed of about 450 km s" " MHansen fc PhinnevI ljl997h con- 
sidered the selection effects, as a result of the flux lim- 
its, of the pulsar surveys and the accuracy of the proper 
motion determinations and found a mean birth speed 
of around 250 — 300kms~^. This is consistent with a 
Maxwellian distribution with (Jk = 190 km s"'^ which has 
a mean velocity of {vk) ~ 30 3 kms ~^. More recently 
lArzoumanian. Chernoff fc Corded (|2002l ) found a best fit- 
ting distribution with two Maxwellian components, one for 



40 per cent of the pulsars with (Tki 
other with (7k2 = 500kms~^. 
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Subsequentl y Hobbs et al.l l|2005l ') extended the work of 
iLvne fc Lorimeil (Il994 ) with a much larger sample of sin- 
gle radio pulsars and claimed that the kick distribution is 
consistent with a single Maxwellian with (Jk = 265 km s^^ 
without a significant low- velocity component. That such a 
kick distribution cannot reproduce the period-eccentricity 
distribution of Be /X-ray binar i es has alrea dy been noticed 
and d iscussed by IPfahl et al.1 (120021 ') and Ivan den Heuvell 
i2003). They concluded that a bimodal kick distribution 
is needed with a low-velocity component ((Tk < 50kms^^) 
when the supernova occurs in a binary system. The idea 
that there might be two types of collapse and therefor e two 
types o f supernova kick wa s origin ally proposed by iKatj 
(|l975h . iPodsiadlowski et al.l l|2004h also find evidence for 
a two-component kick distribution with l ow- ve locity kicks 
mainly in close binaries. iKramer fc Stairj l|2008l ) give an ex- 
tensive discussion of the second kick in the double pulsar 
J0737-3039 and conclude that it p robably had to be s mall. 
Theoretical studies such as those by Scheck et al.l l|2006l ') and 
iKitaura. Janke fc Hillebrandtl (|200(: ) do not yet throw much 
light on this but it may be that there are more lower mass 
supernovae, with smaller kicks, in binary systems or that the 
binary orbit quenches the hydrodynamic instabilities which 
lead to a very asymmetric explosion. 

In a binary system, if the kick is too strong, the sys- 
tem does not remain bound. Even a small velocity kick can 
lead to a large e ccentricity and inclination bet ween the old 
and new orbits (|Brandt fc Podsiadlowskilll995l ). We do not 
know how the angular momentum of the remnant is also af- 
fected by the supernova kick so we do not know how much 
of the misalignment of the neutron star now is caused by the 
orbital inclination. We can however expect the companion 
to continue spinning aligned with the pre-supernova orbit 
immediately after the explosion. There are several binary 
systems with neutron star companions that are observed to 
be misaligned. 



iBrandt fc Podsiadlowskl (|l995D investigated some of 
the effects of high supernova kick velocities on the orbital 
parameters of post-supernova neutron-star binaries. Here we 
look at a variety of velocity kick distributions and consider 
the implications for the distribution of the inclinations be- 
tween the orbit before and after the kick. After ensuring 
that w e can reproduce the work of lBrandt fc Podsiadlowskil 
(| 19951 ) we model Be-star systems with our preferred distribu- 
tions for their progenitors including a range of masses. There 
is somewhat more data now available for comparison in the 
period-eccentricity plane and we find that, though all sys- 
tems can be formed with our models, they tend to be more 
circular than expec ted with kicks distributed according to 
iHobbs et al.l l|2005l ). We investigate what kick distribution 
could lead to the observed eccentricity distribution and also 
how the eccentricity distribution might have changed since 
the supernova. We consider three types of kick distribution, 
a single peaked Maxwellian velocity kick distribution, kicks 
which are direction limited and a double Maxwellian distri- 
bution. 
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Figure 1. The system before the supernova. The two stars of 
mass Ml and M2 are in a circular orbit about their centre of 
mass at the origin. The orbital angular momentum is in the z- 
direction and the orbit is in the xy-plane. When star 2 explodes as 
a supernova it is travelling in the direction of the negative x-axis 
with speed i;orb relative to star 1. It receives a kick of velocity vy^ 
at an angle <j> to the plane of the binary orbit. The angle between 
the projection of the velocity kick on to the binary plane and the 
x-direction is u). 



2 MISALIGNMENT PROBABILITY 
DISTRIBUTION 

In this section we consider the effect of a velocity kick on 
the orbital inclination. We start with a binary in a circular 
orbit. One star then has an asymmetric supernova explosion 
which gives it a kick with velocity ^ Uk < 00 in a direction 
given by the angle (j) out of the binary plane (— 7r/2 ^ (p ^ 
n/2) and an angle between the direction opposite to the 
instantaneous velocity of the star and the projection of the 
velocity kick into the binary orbital plane of ^ a; < 27r 
(see Fig. [1]). For now we assume that no mass is lost. 

We are interested in the misalignment angle of the sys- 
tem, i, after the supernova kick. This is the angle between 
the old and new angular momenta of the orbits (Fig. [21). If 
^ i < 7r/2 then the system is closer to alignment than 
counter-ali gnment and if 7r/2 < t ^ t t it is c loser to counter- 
alignment. [Brand£|&^£odsiadlwsk| l|l995l ) find this angle 
to be given by 



Voih — I'll cos UJ cos ( 



[V^ Sin^ (f) + (f orb — «k COS LJ ( 



(1) 



where iiorb is the initial orbital velocity of the system. The 
relative velocity of the stars after the supernova is 

'^n = I'k + ^orb ~ 2UorbWk COS UJ COS (j}. (2) 

We can rearrange equation ([T]) to find 



Horb 1 I tan( 

COSLJ — ± 



Ilk cos ( 



tan i 



(3) 



when «k 7^ 0, tani 7^ (so that i 7^ 0, tt) and cos(^ 7^ 
(so that (j) 7^ — 7r/2,7r/2). Then if ^ i < 7r/2 from equa- 
tion ([1]) we have «orb > Vk cos uj cos (p. With equation ([Sjl this 



Figure 2. The binary system before and after the supernova in 
the frame of the star 1 of mass Afi. The mass of star 2 is re- 
duced from M2 to in the supernova and it experiences the 
kick. The straight solid arrow is the orbital angular momentum 
of the system before the supernova. The other solid straight line 
is the pre-supernova diameter of the circular orbit indicated by 
the dot-dashed ellipse. The dashed arrow is the orbital angular 
momentum after the supernova. The angle i is the angle of mis- 
alignment between the pre- and post-supernova orbits. If star 1 
has its spin aligned with the pre-supernova orbit, then the angle 
i also measures the misalignment between the post-supernova ec- 
centric orbit and the spin of star 1. The other dashed line is the 
major axis of the post-supernova orbit indicated by a dot-dash 
line. 



corresponds to 



tanj 

and because tani > we see 



tan (h , 
=F r^k cos > 



=F sin > 0. 



(4) 



(5) 



Similarly when 7r/2 < i ^ tt, so that tani < 0, we find 
the same condition as above on sini;/). Now we can rewrite 
equation as 



Vorb 



Wk cos ( 



tan (p\ 
tani 



(6) 



We consider where this equation has real valued solutions 
in the (/) — Hk plane. In Fig. [3] we plot the locus of cos uj = 1 
(^ = 0), 



v+ 



Voih 



1 + 



tan ( 



cos<^ \ tanj 
and the locus of cosoj — ~1 {lj — n), 



"orb 



-1 + 



tan (p 
tani 



(7) 



(8) 



for four values of i. The region between the v+ and «_ con- 
tours is the region where we have real values of cosw and 
it represents the combinations of kick parameters which can 
lead to a misalignment of the chosen i. 
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Figure 3. Possible combinations of the ratio of star 2's supernova kick velocity, tJjj, to the relative orbital velocity, v^rh and angle <j) 
between the kick direction and the orbital plane for four different inclinations i between the pre- and post-supernova orbital planes, top 
left i = 17.2°, top right i = 57.3°, bottom left i = 114.6° and bottom right i = 171.9°. For angles i > 90° the post-supernova orbit 
counter rotates with respect to the spin of star 1. For each value of ip there is a range of values of V]^ which can give rise to the required 
misalignment i. This depends on the angle a) (0 ^ oj < 2it) shown in Fig. [T] In each panel, the solid line corresponds to = (and thus 
to the kick velocity equation [Tjl and the dashed line to ui = n (kick velocity V—, equation [8}. Note that ui = 0, (p = corresponds 
to a kick directly opposed to the motion of star 2 and ui = n, </) = corresponds to a kick in the direction of motion of star 2. In the 
bottom two plots, which correspond to post-supernova retrograde motion, it is not possible to achieve this with a prograde kick (u = it) 
for any value of (p so that the dashed lines corresponding to a; = tt are absent. The dotted lines show the maximum velocity kick, fboundi 
(equation I26I I as a function of (p, for which the system remains bound. Below this line the system remains bound after the supernova 
and above it the binary is disrupted. In order to find the probability distribution P{i) of the misalignment angle i in equation l|18|l we 
integrate in (Di^/tJorb i 0)-space over the regions R = Ri + R2 between the contours of (solid line), V— (dashed line). To find only the 
probability distribution P{i) for the bound systems alone we integrate over the regions R2 only. 



3 ISOTROPIC MAXWELLIAN KICK tributed over a sphere so that the angles defined in Fig. \T\ 

DISTRIBUTION are distributed as 

As an illustration we apply these results to a simple isotropic 

Maxwellian kick distribution. For an isotropic kick distri- ^ 

bution the direction of the kick velocity is uniformly dis- P{Lo)du — — dco (9) 
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and 



= cos (f> d(j). 



(10) 



We here choose the kick speed to have a Maxwelhan distri- 
bution so that 



P(iik)duk = 



2vi 



(11) 



where ak is the dispersion of the velocity, and recaU that 
iHobbs et all (2005) find (Jk = 265kms"\ In sections K2 
and 15.31 we shaU consider ahernative velocity kick distribu- 
tions. Because i = i(v\^,(f>,uj) its probability distribution is 



P{i)di = 



fk=0 J 4>= 



P{(l))P{uj)P{vi,) duj d4>dv^ 



(12) 



We change variables from {v-^^ 4>^^) to (fk, i) and find 

J J P{(l>)P{uj)P{vy)\J\d4>dV],di, (13) 

where R = R\ + R2 is, the region in the {(j), Wk) plane where 
cosoj is real valued. This is illustrated in Fig. [3] for different 
values of i. Outside of the region bounded by these curves, 
the given velocity kick and angle cannot produce a system 
misaligned by i because then | costj| > 1. We consider this 
region for bound systems in Section 13.11 The Jacobian, J, 
for the change of variables is given by 



duk d(j) duj = \J\ d«k d(j> di, 



where 



J : 



9vk 

dvy. 



8vk 



\9<t> I 
\84> 1 



cK 
(t) 
(t) 



1'-k,<t> 



(14) 



(15) 



Because ui, (j) ^nd ?;k are independently distributed we find 

(16) 



J = 



'di 



"k.'^> 



We differentiate equation (|6]) to find 

' dio \ I tan 

Yi 



J 



sm Lu sm t 



(17) 



Now we can write down the full probability density distri- 
bution for the misalignment as an integral 



P{i)di 



2al 



idwk di. 



where 



2 

WkB 



svaLo \ sm i 



(18) 



(19) 



and uj — uj{i,V]^) is defined by equation ((6)1. 

This can now be integrated numerically for specific val- 
ues of Vorb/o'k to find the probability distribution of the 
misalignment angle, i, of the system. We plot this as the 
solid lines in Fig. |3] for Horb = fk (upper line on y-axis) 
and Horb ~ 0.5 (Tk. The higher the misalignment angle, i, 
the lower its probability of forming. The higher iiorb/fk, the 
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Figure 4. The probability distribution P{i) of the post- 
supernova misalignment angle i is plotted as a function of i for 
the case when the probability distribution of the kick velocity 
Vk is a Maxwellian with dispersion crk and no mass is lost from 
the system (/ = 0). The solid lines correspond to integrating 
over the full area R = Ri + R2 in Fig. |3] and so to all systems 
whether or not they remain bound. The curve with the higher 
value of P(i = 0) corresponds to crk/tiorb = Ij and the other one 
to o'k/^'orb = 2. These are normalised so that f^^^ P(i)di = 1. 
Note that the distribution with the larger velocity spread is more 
able to produce counter-rotating systems. The dashed lines are 
the corresponding probability distributions (ck/^orb = li upper; 
(Jk/forb = 2, lower) with the integral now taken over only the 
region R2 in Fig. [3] and so only including those systems which 
remain bound. The change in P{i) for given values of i and (Tk is 
due to the fraction of systems which become unbound. The dot- 
ted lines are the probability distributions P(i) for systems which 
remain bound, with Ck/norb = li but for which mass loss also 
occurs. The fractional mass loss is / = 0.2 for the upper dotted 
line and / = 0.5 for the lower. 



lower the relative probability of a counter-rotating system 
and the higher the probability of a system being close to 
alignment. 



3.1 Bound Systems 

We have calculated the probability distribution of the mis- 
alignment angle of the system but have not yet considered 
whether the resulting system remains bound. If the kinetic 
energy of the new system exceeds the gravitational poten- 
tial energy then the system does not remain bound. We no 
longer have a binary system and the two stars fiy apart. 
Initially the gravitational energy of the system is 

GM1M2 



and the kinetic energy is 



1 /Ml Ma 



2 

^orb 5 



(20) 



(21) 
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in the centre of mass frame, where M = Mi + Al2. In a 
circular orbit we have 



and the total energy 



Eto 



1 GM1M2 



2 a 



So 



_ _ 2 Ml Ma 

-Cigrav — forb 



(22) 
(23) 
(24) 



After the kick the gravitational energy remains the same 
because in this case no mass is lost. The kinetic energy be- 
comes 



, ^1 / Ml Ma ^ a 



(25) 



where Vn is given by equation ([2}. If -E^in > — -Egrav the new 
system is unbound. The condition for the system to be un- 
bound is > 2v^,.^^. We s( 
for the critical velocity of 



bound is > 2v^,.^^. We solve Vn = 2v^^y, with equation ([2]) 



Abound — 



«orb sm ( 



tanz 



± 



3«orb +«orb^ 



(26) 



Because Hk > we take the term with the positive sign. If 
t^k > Abound then the system is unbound but if «k < Wbound 
it remains bound after the supernova kick. 

This condition for the system to be bound affects the 
region in the Wk-0 plane that we integrate over to find the 
probability distribution for the misalignment angle. To in- 
tegrate over all systems, bound or unbound, we integrated 
equation (|18p over the region R — Ri + R2 shown in Fig. |3l 
There we also plot the upper limits on the velocity kicks for 
the system to remain bound as dotted lines, Hbound- Below 
these lines a system remains bound but above the kick is 
too strong and the two stars fly apart. The region in which 
we have bound systems, R2, is much smaller than the region 
that can produce the given misalignment angle, _Ri -I- i?a . To 
find the probability distribution of the misalignment angle, 
P{i), for bound systems only we integrate expression psp 
over this smaller region 7? = _Ra . 

We compute this numerically and plot it in Fig. |4] as 
the upper dashed line when iiorb = fk- The lower dashed 
line is for Uorb = O.Sffk- Most kicks unbind the systems but 
those that remain bound are somewhat more likely to be 
counter-aligned than for smaller kicks. On the other hand if 
Voih > o"k we find that few kicks are able to cause counter- 
alignment. 

As expected the higher the misalignment angle of a sys- 
tem the lower the probability of it forming. We see that, by 
restricting to only bound systems, the number with small 
misalignment is greatly reduced whereas those closer to 
counter- alignment are less so. The probability of a system 
closer to counter-alignment than alignment becomes rela- 
tively high for small Voib- 

3.2 Mass Loss 

So far we have assumed that the mass lost from the system 
is negligible. We now allow the mass of star 2 to fall in the 
supernova to 



Ma = Ma - /M, 



(27) 



where / is the fraction of mass lost relative to the total mass 
of the binary system so that the total mass of the system 
becomes 



M' = (1 - f)M. 



(28) 



The gravitational energy of the system after the supernova 
is 



'a Ml Ma 



rpl _ _ '2 



(29) 



where v'^^.^, is the relative velocity that the stars would have 
in a circular orbit of separation a and 



«orb = 



= ^orb-^ = (1 - /)t'orb 



(30) 



because the instantaneous separation a' — a. Note that v'^^^, 
is the orbital velocity for the equivalent circular orbit while 
the orbit itself is eccentric after the supernova. The new 
kinetic energy is 



, ^ 1 2MjM 
-e-kin M' 



(31) 



and so the condition for a bound system (-B^in < ~£'grav) 
becomes 



"n < 2(1 - /)Uorb- 



Thus we find 



Abound 



Horbl sin(j!)| 



(32) 



(33) 



tani ■ V ■ """^nan^r 

In Fig. [5] we plot contours of Abound for varying / with i — 
0.3 rad = 17°. 2 and Vorh = c"k- The and w_ contours and 
the top dotted line remain the same as in the top left plot 
in Fig. [3] The more mass that is lost in the supernova, the 
lower is the limit on the kick velocity for a bound system 
and so the less likely a bound system with a given inclination 
becomes. 

In Fig.|3]for ~ we plot the probability distribution 
for bound systems for f — (upper dashed line), 0.2 (upper 
dotted line) and 0.5 (lower dotted line). A larger / increases 
the likelihood of counteralignment in bound systems. 



4 ECCENTRICITY PROBABILITY 
DISTRIBUTION 

In the previous section we discussed the effect of particular 
supernova kick distributions on the distribution of orbital 
misalignments. We now consider what kick distributions are 
most able to give rise to the observed Be-star eccentricity 
distribution. Given the sparsity of the data, the large num- 
ber of free parameters and the unknown selection effects, we 
do not attempt to find a best fit to the periods and eccen- 
tricities of Be stars. Rather we look for a kick distribution 
consistent with these observations and then examine its con- 
sequences for the distribution of orbital misalignments. 

The new semi-major axis of the orbit after the super- 
nova can be found from 



a an 



(34) 



where a, the old semi-major axis, is the instantaneous sep- 
aration. Combining this with equation ((2} we can find an. 
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Figure 5. This corresponds to the top left panel in Fig. |3] for 
i = 17.2°. The solid line and the dashed line are the same in 
both Figures. The dotted line in Fig.|3]is the same as the upper- 
most dotted line here and represents the upper boundary for vy^ 
which permits the post-supernova system to remain bound when 
no mass is lost. When mass is lost from the system the upper 
limit on V]^ such that the system remains bound decreases as the 
fraction of mass lost / increases. The dotted lines here correspond 
to / = 0, 0.2, 0.5 and 0.8. 



The new system has specific angular momentum 

h' = r X Vn. (35) 
where r is the separation vector of the stars. We have 



GAf'a„(l - e^) = \r x 



(36) 



and so 
GM'an(l 



r 2 . 2 

V]f^ sm ( 



' -I- (t;k cos LO cos ( 



■ Uorb)^] 



(37) 



which c an be solved to find the eccen tricity, e, of the new 
system (|Brandt fc Podsiadlowskilll995i ). The binary system 
is unbound if e > 1, in which case an ^ 0. Although we 
could find eccentricity probabilities by direct integration in 
a similar way to the inclinations in the previous section, it 
becomes very complicated and we do not learn much new 
from the procedure. Instead we use Monte-Carlo methods 
to evaluate the integrals. 

We note that e and i both depend only on Wk/worb, 
1 — / and the two angles and u. For typical progenitors 
of Be stars the dependence on masses, through Uorb and 
1 — / as well as the radius of star 2, turns out to be weak 
so it can suffice to concentrate on only one set of masses 
initially. We choose a pre-supernova mass of M2 = 5 A4q 
that leaves a neutron star of mass M2 = 1.4 Mq and a 
companion mass Mi = ISMp. These masses were used by 
iBrandt &: Podsiadlowskil (|l995l ) and we have ensured that 
we can reproduce their results too. 

We use the NAG Library routine G05CAF to generate 



pseudorandom numbers {Xi } uniformly distributed between 
and 1. Then for an isotropic kick distribution 



Xji 



and 



J2- 



(38) 



(39) 



The distribution of orbital periods. Pi, immediately be- 
fore the supernova depends in a complex way on the pre- 
vious evolution of the system. There are many as yet 
unquantified processes tha t contribute to this evolution 
l|Hurlev. Tout fc Pols|[20o3 ) and so we stick with the rela- 
tively simple assumption that log Pi is uniformly distributed 
between Pmin and Pmax so that 



l0gPi = l0gPi-f X,3(l0gPn: 



log Pmin 



(40) 



We take Pmin to be the period at which star 1 would fill its 
Roche lobe in a circular orbit if it has the main-sequence 
radius given by I Tout et al.l ([l996', 5 Rq for a 15 Mq star) . 
Its Roche - lobe radius J?l is approximated by the formula of 
lEggletonl digsl). 



Pl 



0.49(7! 



0.6(73 + logg(l +qc 



< (7 < 00, (41) 



where q = M1/M2. We take Pmax = 10^ d because beyond 
this almost all systems are disrupted. 

For the post- supernova systems 

iBrandt fc Podsiadlowskil (|l995t ) rejected any system 
that would have filled its Roche lobe if it were circular at 
its periastron separation. So if 



(1 



. Pi 

oil) 



(42) 



they rejected the system. In practice we expect that sys- 
tems cannot actually survive down to this separation be- 
cause ti des enforc e pseudosynchronization of star 1 at peri- 
astron (lHutlll99ll) and so it ends up spinning up to about 
1.16 times faster than it would in a circular orbit of the peri- 
astron separation. However there is no equivalent potential 
theory in the eccentric orbit so we do not try to be any more 
precise than condition (|42p . 

We can reprod uce figs 4 and 5 of 

IBrandt fc Podsiadlowskl (|l995l ). They chose «k to be 
constant and used a period distribution which is uniform 
in Pi rather than logP. In their figs 4-6, for a given a;-axis 
value, they found the median on the y-axis of 10,000 runs 
and then worked out the regions in which 20, 40, 60, 80 
and 98% of systems lie away from that median. They found 
very high values of the mean inclination because they used 
a high single value kick velocity. 

Instead we integrate over Wk, distributed according to 
equation (|lip up to 850kms~^ using Simpson's rule and 
ui, cj> and P by the Monte Carlo method. In our figures we 
prefer to plot contours of probability density in the 2D space 
normalised so that the probability of lying in the plots is 1. 

Rather than sticking to the fixed masses we distribute 
the masses of the companion star from Mimin = 5 M© to 
Afimax = 25 Mq accordiug to a mass function 

N[Mi) dMi oc Mi"^-'' dMi (43) 

(|Kroupa. Tout fc Gilinor3 1 19931 ) which can be generated 
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from 

M. = (^i^)"'^, (44) 

where we find Xq and k from the minimum and maxi- 
mum masses. Because the mass range is hmited the re- 
suhing distributions are not very different from the fixed 
ini tial masses of Mi = 15 an d M2 = 5Mq as used 
by iBrandt fc Podsiadlowskil (|l995l ) to represent a typical 
Be star binary. Before presenting our results we discuss the 
observations with which we compare. 



5 MODEL COMPARISON TO OBSERVED 
SYSTEMS 

There is now a large number of Be star binary systems with 
measured periods and eccentricities. The data to which we 
shall apply our models are given in Table[T] The bulk of these 
come from the catalogue o f Be/X-ray binaries assembled by 
iRaguzova fc Popovl l|2005l 'n. We also include O /B stars that 
we have found in the literature that have no emission but 
must have formed in the same way. We can use these to look 
at the eccentricity distribution and remember that the disc 
in Be stars may come and go so that B stars can become Be 
stars and vice versa. 

Because we are primarily interested in Be stars with 
discs we also impose the condition that there must be enough 
room for a decretion disc of 4 R* or so inside the Roche lobe 
of star 1. We explain our choice in the appendix but note 
that there is one Be star in Table [T] 0535-668, that we do 
not include in our analysis because it cannot accommodate 
such a disc. 

There are three systems for which the misalignment an- 
gle i between the Be star spin and disc have been estimated. 
There is one more system, 28 Tau, which has no reliable or- 
bital parameters but for which an inclination of i = 25° —30° 
has been suggested. 

5.1 Single Maxwellian Peak Velocity Kick 
Distribution 

We now co nsider how well the standard velocity kick dis- 
tribution of iHobbs et all (|2005l ') fits the observed data. In 
the left frame of Fig. [6] we plot the probability distribu- 
tion contours of eccentricity against final period for these 
systems and in the right panel we plot contours in the 
inclination-eccentricity plane. The probability of a highly in- 
clined system is quite small and most systems end up closer 
to alignment than to counter-alignment. Larger inclinations 
are more likely for the more eccentric systems. 

To illustrate this further, the solid lines in Fig. [7] show 
the probabilities of e and i integrated over all systems. We 
also plot a histogram of the observed Be star systems' ec- 
centricities. It is evident, in line w ith previous findi ngs on 
less substantial data sets, that the iHobbs et al. I 12005) dis- 
tribution gives a poor fit to the eccentricity distribution. 
It does, however, give a fairly fiat distribution of misalign- 
ments. There are too few measured misalignments to use 

^ http:/ /xray.sai. msu.ru/~raguzova/BeXcat 




e 



Figure 8. Solid lines are the cumulative eccentricity probability 
distribution for (Tj; = 265, 15, and 5kms~^ in order of increasing 
height on the right hand side in the plot. The cumulative eccen- 
tricity probability distribution for the observed systems is shown 
for the twenty Be stars (dotted line) , the ten B stars (dashed line) 
and for the total of thirty B and Be stars (solid line). We use a KS 
test to compare the observed distribution of system eccentricities 
to our model predictions with the results given in Table [2] The 
quality of fit peaks for the Be stars around (T]j = 5kms~'^ and for 
the combined Be and B stars at around (Tj^ = 15kms~^. 

these as a test b ut the measu r ed va lues of are consistent 
with the standard iHobbs et al.l (|2005l ) kick distribution. 

Because the final period distribution depends very 
strongly on our choice of initial period we have not made 
use of it for a statistical comparison with observations ex- 
cept to say that the systems appear to fit well in Fig. [S] 
(and later in Fig. llOf) . This implies that our choice of period 
distribution is reasonable. 

We use the Kolmogorov-Smirnov test (KS test) to de- 
termine if two datasets are significantly different. It is non- 
parametric and distribution free. In Fig. [8] we plot the cumu- 
lative eccentricity distributions of the Be stars, of the B stars 
and of the B and Be stars combined. We also plot cumula- 
tive eccentricity distribution for our model predictions. We 
then perform a KS test between the observed data and our 
model prediction and give the results in Tabled We find the 
largest deviation of the observed data from the model and 
use probability tables for the KS test to find the probability 
that the observed sample of stars came from the distribution 
predicted by the model. In Table [2] we give the probability 
that the observations of Be stars. Pi, B stars, P2, and com- 
bined B and Be stars, P3 are consistent with our various 
models. Probabilities smaller than 10"'* are li sted as zero. 
The very small probabilities for the standard iHobbs et al.l 
(2005) distribution (the first line in Table [2| demonstrates 
that it is essentially impossible that the B and Be systems 
have formed as they are in this way, with a single Maxwellian 
distribution with ak — 265 km s"'^. 

We cannot, however, immediately rule out a kick distri- 
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log,„ (P,/d) ^ 

Figure 6. Left: the stars listed in Table [T] in a period-eccentricity diagram. The solid squares are the Be/X-ray binaries and the open 
triangles are the binary B stars. The contours are lines of constant probability density P(Pf,e) when the dispersion of the velocity kick 
distribution in O]^ = 265 km s~^ and post-supernova binaries which are too tight to permit a disc of size 4_R* are excluded (section |4]l . 
The probability density P is defined so that the probability of finding a system with e in the interval (e, e + de) and with period log Pf 
in the interval (logPf,logPf + dlogPf) is PdedlogPf. The area at small Pf is excluded by the models. The outermost contour is at 
P = 0.01. Moving inwards the contour levels are P = 0.1, 0.2, 0.4, 0.6, 0.8 and 1. Right: contours of equal probability density P(e,i) 
for the same models as in the left panel. The contour levels are, starting from the left, P = 0.01, 0.1, 0.2, 0.4 and 0.6. The asterisks 
correspond to the systems given in Table [T] for which there are estimates of the misalignment angle i. 




Figure 7. Left: the two dimensional probability density P{P[,e) contours, of which are shown in Fig. [6] (left), integrated over period 
to give a one dimensional probability density P(e). This is plotted as the solid curve (normalised so that Jq P(e) de = 1) and for a kick 
distribution with (Tk = 265 km s~^. The dotted line gives the probability density P(e) for models constructed with a-]j = 15 km s~^. The 
histogram is the observed eccentricity distribution for the Be stars (the solid squares in Figure [gjl. The lower value of cr]j gives a much 
better fit to the data. Right: as for the left panel but the probability distributions for misalignment angles i predicted by the models 
are plotted. The few measured misalignment angles which are known are not plotted but all exceed i ^ 25° . It is evident that the solid 
line ((Tj; = 265 km s~^) is consistent with a fairly uniform spread of misalignments while the models with ctjj = 15 km s~^ (dotted line), 
which provide a better fit to the eccentricity distribution, are not consistent with the observed misalignments. For this model only 7 per 
cent have i > 15° and only 1 per cent have i > 25° . 
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Spectral Type 


Pf/d 


c 


i 


Systems with Emission 


0053+604 (7 Cas) 


B0.5 IVe 


203.59 


0.261 


about 25° 43 


0115+634 


B0.2 Ve 


24.3 


0.342 




0331+530 (BQ Cam) 


08-9 Ve 


34.3 


0.33 




0352+309 (X-Per) 


09.5 IIIe-BO Ve 


250 


0.11* 




0535+262 (V725 Tau) 


BO Ill-Ve 


111 


0.475 




0834-430 


BO-2 Ill-Ve 


105.8 


0.12<5 




J1008-57 


09e-Ble38 


247.5 


0.66"'° 




1417-624 


Bl Ve 


42.12 


0.446^ 




1845-024 


Be 


242.18 


0.88^'' 




J1946+274 


BO-1 IV- Ve 


169.2 


0.33" 




J1948+32 


BO Ve 


40.415 + 0.010 


0.033 + 0.0131'! 




2030+375 


BOe 


46.0202 


0.41610 




J2103. 5+4545 


BO Ve 


12.66536+0.00088 


0.4055+0.0032" 




SAX J0635.2+0533 


B2V-BlIIIe 


11.2 ± 0.5 


0.29 + 0.0912 




XTE J0421+560 


B4 III-V[e] 


19.410 


0.6213 




4U 2206+543 


09.5Ve 


9.570 


O.I515 




B1259-63 


Be 


1236.72404 


0.86988691*^ 


greater than 55° 42.16 


* 0535-668 


B2 Ille^S 


16.65 


0.82+0.041'' 




IE 1145.1-6141 


B2Iae 


14.365 


0.20^0 




GRO J1750-27 


Be 


29.817 


0.360^6 




LS I +61 303 


BO Ve 


26.5 


0.63±0.11''i 




Systems without Emission 


1901+03 


OB 


22.58 


0.036'* 




J0045-7319 


B 


51.16926 


0.807981** 


25 - 41° 16 


J1740-3052 


B 19 


231.02965 


0.5788720 




2S 0114+650 (LS I +65 010) 


BO. 5 lb 


11.600 


0.18 + 0.0521 




* 4U 1538-52 


BO lab 


3.730 


0.18 + 0.0122 




4U 1907+09 


08-9 la 


8.380 


0.22+0.0523 




BP Cru (GX 301-2) 


B1.5Ia 


41.5 


0.4624 




GAG 1657-415 


BO-6 lab 


10.44809^5 


0.10426 




Vela X-1 


B0.5 lb 


8.964368 


0.08982'' 




* LMC X-4 


08 III 


1.40841 


0.0062* 




J1638-4725 




1940.9 


0.95529 




XTE J1855-026 


or B 


6.067 


0.04 + 0.0232 





Table 1. Binary Be stars and binary B stars that we use in our analysis. We include all the binary B and Be stars we have been 
able to find reference to in the literature for which both period s and ecce ntricities are re liably known. Most of the data that we 
use comes from the catalogue of Be/X-ray binaries assembled by iRaguzova & PopoyI l l2005h . We also include measurements or esti- 
mates of the misalignment angle i for the few systems for which it is available. The systems marked with an asterisk (*) are not 
included in our analysis because th e stellar separatio n is not large e nough to accommod ate a B e -type disc of size 4 . Refe rences 
for quantities in the Tab le are given l[Harmanec et al.l Jio OOl ^Rappapo rt et alj 17978') 3stcll^^^J^98^}jDelgado-Martf ct al l 1I2OOII') 

J Negueruela et al.l 1I2OOOI') IWilson et al.Nl997l') iFinger.^Vilson fc Ghakrabartvl lll99a) teallowav. Wang fc Morgan (2005) " Wilson et al. 



I2OO3I') 1*1 Wilson. Fabregat fc Cobu rn (2005) "iBavkal et al.l ||2007|) l^Kaaret. Cusumano fc Saccd [2 0001 i-' Barsukova et al.l 
" Gallowav. Morgan fc Levinj i2004) i ^Ribo et al. | ||2006'1 I'^Hu ghes fc Bailed lll999l) llHutchings et a l. (1985) i*^ Kaspi et al.) 

5rS 



Stairs et 



12001 



'iRay fc Chakra b^tvl [2002 ) 2^Grundstr om et al.l 72007. 

van Gendercn & Stcrkcn' ^2007) 25Bildsten ct a.1: H sOTt) 2'^Cha.k rabartv et al 
McLaughlin (2004) 32 Corbet & Mukai (2002)_ 35 Ch arles et al.l ((19 83) 36 Scott et 
Gkazaki fc Negueruelal ||200l|) 4llC as ares 



^2006) 



^ M akishima et al.l 



iRaichur et al.l 
(I1993I') 2^ ildsten c TaLl (|1997|) 2fLevine et al.[ 

. .^ _ et al.l il9m ) 37 Fing er et al.l l ll999l) 3 ^Coe et al 

et al.l l(2005h 4 ^Wex et al.l [l99^ and 4i iHummell ((19981) . 



2007 
19941') 
19841) 
1991) 
I994I ) 



bution in the form of a single Maxwellian distribution with 
aic = 265 km s"'^, because the current eccentricity distribu- 
tion might not be representative of the eccentricity distri- 
bution immediately after the supernova. One possibility is 
that the systems have begun to circularise by some mech- 
anism but, given that they are not completely circular, the 
timescale on which this circularisation operates is coinci- 
dentally close to the time since the supernova. We return to 
this in Section [S] A second possibility is that there might be 
a selection effect on the observed systems. The interaction 
that leads to emission may only occur close to periastron 



in the very eccentric systems that spend very little of their 
time there. Alternatively it may just be that periods and ec- 
centricities are more easily measured for the least eccentric 
systems. We again return to this in Section (6) 

On the other hand the observed eccentricity distribution 
is consistent with a Maxwellian kick distribution provided 
that the kick distribution is peaked at a lower velocity. In 
Table[2]we show the effect of reducing ak and find very good 
fits to all observations when 10 < cr^ /km s~^ < 20. In Fig. [7] 
we also plot the dotted lines to show the eccentricity and 
inclination distributions for a single Maxwellian peak with 
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Single Maxwellian Distribution 



and 





Pi 


P2 


P-i 






265 


0.0003 


0.0025 


0.0000 






190 


0.0003 


0.0025 


0.0000 






80 


0.0013 


0.0035 


0.0001 






40 


0.0199 


0.0088 


0.0057 






20 


0.1572 


0.0217 


0.1000 






15 


0.2550 


0.0238 


0.1308 






10 


0.3840 


0.0213 


0.1141 






5 


0.4295 


0.0103 


0.0445 






2 


0.3699 


0.0091 


0.0378 






Bimodal Maxwellian Distribution 








0"k2 


Wl 


Pi 


P2 


Pa 


15 


265 


0.4 


0.2087 


0.0217 


0.1168 


15 


190 


0.5 


0.0210 


0.0088 


0.0061 


90 


500 


0.4 


0.0009 


0.0032 


0.0000 


15 


500 


0.4 


0.2475 


0.0229 


0.1279 



Table 2. The KS test probability values, the likelihood that the 
observed systems are chosen from the model distribution. In the 
upper part of the Table the velocity kick distribution is modelled 
as a single Maxwellian. In the lower part of the Table the velocity 
kick distribution is modelled as two Maxwellians with o-jji and 
CTii2 and Wl the relative weighting of that with a]g^i to that with 
CTk2- The likelihood for the Be stars only is Pi, B stars is P2 and 
for the combined B and Be stars is P3. Probabilities less than 



10" 



appear as zero. 



dispersion (Tk = 15kms~^. We see that this curve appears to 
fit the eccentricity distribution much better than that with 
ak = 265kms~^. However, for such low velocity kicks, the 
misalignments tend to be small (right panel of Fig. O and 
are hard to reconcile with the observed values. 



5.2 On Axis Kicks 

It is possible that the direction of the velocity kick in the 
supernova is restricted l|Brandt fc Podsiadlowskilll995l ). To 
investigate this we consider the extreme case that the kick is 
always directed along the spin axis of the star, perpendicular 
to the orbital plane, so that (p = 7r/2 and uj is undetermined. 
With equation (|28|) and v'^^.^, = GM/a, we can express equa- 
tion (l34| as 



vl = Uorb(l - /) 2 - 



and then equation (|37l) becomes 



(45) 



^^k • 2 
-T5— Sin 

orb 



■ cos LU COS ( 



Worb 



.(46) 



If the kick lies on the z-axis we have sin (j) — 1 and so the 
misalignment angle between the old and new orbital planes 
is 

1 



(47) 



+ 1 



2 2 2 

"11 = +Voih- 



(48) 



We can then relate the eccentricity and the inclination by 
2 1 



1-/ 



2 - 



/ 



(49) 



and we plot this in the right panel of Fig. |9] For A42 — 5 Mq 
and M2 = 1.4 Mq we have 



■' ^ M 



1.4 + Ml /Me 



(50) 



5 + M1/M0 

For illustration we choose Mi = 5, 15 and 25 Mq and so 
/ = 0.36, 0.18 and 0.12. 

We plot the eccentricity distributions if the kick is par- 
allel to the binary orbital axis in the left panel of Fig. |§] We 
see that we cannot get highly misaligned systems or low- 
eccentricity systems with a kick in the z-direction. The low 
eccentricities could be explained by circularisation but the 
highly misaligned systems are ruled out in this case. 

5.3 Bimodal Distribution of Velocities 

A bimodal velocity kick distribution has been suggested to 
explain both the high- velocity neutron stars and also the fact 
that neutron stars appear to be easily contained in globular 
clusters (Katz 1975). The escape velocity of a neutron star 
from a globular cluster is about 30 kms~^ and it is generally 
believed that about 10% of neutron stars born within them 

are retained _^rukicr l oji). 

lArzoumalian. Cher'^ff fc Corded (l20oi ) used observed 
properties of radio pulsars and other neutron stars to show 
that a two component velocity distribution fits the data 
much better than any one component model. They used a 
distribution of velocities with two Maxwellian distributions 



P{vk) dvk (xwi 



+ {l-wi)^e 



' dvk 



(51) 



where Oki ~ 90kms~^ and ak2 ~ 500 km s^^ and wi — 
0.4 ± 0.2 is the weight of the first distribution. We perform 
the same KS test on this bimodal distribution and find that 
it does not fit our eccentricity distribution well at all (see 
Table [2]). We note that almost all systems that fall into the 
higher peak are disrupted so that the poor fit is entirely due 
to the high dki- 

We found previously that the best fitting distribution 
with one peak had ak — 15kms~^. For a bimodal distribu- 
tion with Oki = 15kms~^ and crk2 = 265kms~^ with equal 
weight we find that the fit is somewhat poorer than for the 
single peak at Ok ~ 15 km s~^ . However it fits the data signif- 
icantly better than a distribution with just ak = 265 kms~^. 

The weight factor wi is unimportant if (Jk2 is large 
enough to disrupt most systems. For aki — 15kms~^ and 
(Jk2 ~ 500 km s"'^ with wi — 0.4 the fit is as good 
as the single peak at Ok = 15kms~^ (see Table [2]). In 
Fig. [TO]we plot the eccentricity-final period and inclination- 
eccentricity contours for this velocity kick distribution. Most 
of the surviving systems are from the o"ki = 15kms~^ part 
and so the inclinations are all low. Nearly all systems with 
kicks from CTk2 ~ 500kms~^ are disrupted. 
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log,„ (P,/d) 



Figure 9. Left: as for Figure [6] except tliat tlie direction of the velocity kicks are confined to the 2-axis (the direction of orbital angular 
momentum in the pre-supernova binary). As in Figure|S]the contours of probability density range from P = 0.01 to P = 1. The solid 
squares are the Be stars and the triangles are the B stars without emission from Table [T] The models cannot account for the observed 
low values of the eccentricity in some systems. Right: in this case, because the kick direction is fixed, there is a simple relation between 
misalignment angle i and eccentricity e for a given fraction / of mass lost in the supernova (equation I49I I. This relation is plotted for 
/ = 0.12 (solid line), / = 0.18 (dotted line) and / = 0.36 (dashed line). The stars show data for the systems which have estimated 
misalignment angles. 



This bimodal distribution fits the distribution of eccen- 
tricities we find here just because the low-kick systems dom- 
inate those that remain bound. The effect of the larger kicks 
is minor simply because most systems with large kicks are 
disrupted but this would explain the high pulsar space veloc- 
ities. Until we understand the nature of the supernovae ex- 
plosions properly in 3-D, bimodal distributions of this kind 
will remain good but ad hoc solutions. 



6 DISCUSSION 

We have modelled the eccentricity distribution of observed 
B and Be stars and we confirm the findings of others that 
the e ccentricities t e nd to be lower than predicted by a stan- 
dard iHobbs et al.l (|2005l ) kick distribution. The final peri- 
ods depend on the uncertain choice of inner boundary and 
initial period distribution. So, to test the models, we com- 
pared only the probability distribution of the eccentricities 
(left panel of Fig. ■ To quantify this we performed a K-S 
test (see Table [2| to see how well the B and Be stars fit 
our models. To obtain a good fit from a Maxwellian kick 
velocity distribution we requi re Jk = 15kms~^, m uch lower 
than the (Jk = 265 km s"^ of iHobbs et~al] (120051 ). This can 
be combined with a second Maxwellian as long as that has 
ffk large enough to disrupt most systems. A combined low 
and high velocity kicks in the bimodal distribution can then 
reproduce the high space velocity pulsars and the eccentric- 
ity distribution of the bound systems. 

In Fig. [7] we show that the few estimates of misalign- 
ment angles i between Be star spin axis and disc axis are 
too large to be easily reconciled with kicks with CTk as 



low as 15kms~^. They are much better accommodated by 
(Jk — 265kms~^ (Fig. [6)|. This is even more true for the B- 
star binary PSR J0045-7139 in which the B star spins ret- 
rogradely with respect to the orbit. However, although the 
observed misalignments are indicative of high kick velocities, 
we cannot yet base any firm conclusions on the inclinations 
because (i) there are only a few, rather uncertain, measure- 
ments, (ii) it is easier to measure a large inclination than 
a small one and (iii) misalignment angles of les s than the 
disc opening angle, estimated to be around 13° l|Hanuschikl 
Il996l ). would not be sufficient to easily give rise to change 
between Be star and shell star. 

If on the other hand we beli eve that Jk must 
be larger, as in the distribution of iHobbs et all (|2005l . 
Uk = 265kms~^) and even the bimodal distribution of 
lArzoumanian. Chernoff fc Corded ()200l CTki = 90 kms , 
f"k2 ~ 500kms~^) then we need an alternative explanation 
of the lower than expected eccentricities. We first consider 
whether this could be due to observational selection effects. 
We are less likely to observe systems with high eccentricity 
if we only see them when they are close to periastron. Such 
systems spend the majority of their orbital periods closer to 
apastron and so are less likely to have been observed. In the 
case of the Be stars this may simply be due to the fact that 
the companion is much more likely to interact with the de- 
cretion disc at periastron. The same would not be true of the 
plain B stars which appear to follow a similar eccentricity 
distribution. However there are fewer of them and they still 
may have been selected by radial velocity variations which 
would be larger at periastron. 

If we assume that we have found all of the Be star sys- 
tems in our lowest eccentricity bin of < e < 0.105 (see 
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Figure 10. As for Fig. [6] except that tlie model velocity kick distribution consists of two Maxwellians with a]^i = 15kms~^ and 
f7i;2 = 500kms~^ and relative fraction of the former wi = 0.4. The contours in both plots are as in Fig. |5]with contours ranging from 
P = 0.01 to P = 1. Left: the contours of constant probability density in the eccentricity-final period plane. Solid squares are the Be stars 
and triangles are the B stars without emission in Table [T] By comparison with Fig. |5]it is evident that this model provides a better fit 
to the observations. This is confirmed by the results of the K-S test given in Table [2] Right: the probability density in the (e, i)-plane. 
The asterisks are the systems which have an observed misalignment angle. The observational points are at low values of the probability 
density, so the observed misalignment angles are hard to reconcile with such a model. 



histogram in Fig. [7]) then the relative number needs to be 
brought down by a factor of about 3.4. This would require 
the actual number of Be star systems in our galaxy to be 
3.4 times greater than what we have found and the extra 
ones must all have high eccentricities. We have found mea- 
sured eccentricities and periods for 20 Be systems in our 
galaxy and so would need 3.4 x 20 = 68 systems to ac- 
count for selection effects. There are 67 observed Be star 
systems in our galaxy in the Be/X-ray binary catalogue 
jRaguzova fc Popovl [2OO5I ) of which 52 do not have both 
a measured eccentricity and period. It is unlikely that the 
Be stars without observed eccentricities will all turn out to 
have high eccentricity. However it is important for observers 
to measure the eccentricities of more Be star systems be- 
cause this is vital to rule out selection effects. 

Alternatively, and perhaps more interestingly, it is pos- 
sible that Be stars did form with an eccentricity distribution 
favouring large eccentricities of e > 0.5 or so. In this case 
the systems must have subsequently circularised and the cir- 
cularisation timescale must be similar to the lifetime of the 
Be phase. The periods of Be stars can be large and the tidal 
circularisation time, for these stars with radiative envelopes, 
is much longer than their lifetimes. One possibility is that 
the neutron star interacts with the decretion disc of the B 
star at periastron passage. Such a mechanism might be self 
regulating in the sense that the decretion disc only has time 
to build up to a large radius when the system is very eccen- 
tric and the neutron star spends a long time far from the 
B star. So this merits further investigation particularly if it 
can be established that supernovae kicks cannot have such 
a low dispersion of (Jk = 15kms~^. 



7 CONCLUSIONS 

The distribution of eccentricities in Be and B stars cannot 
be reproduced directly if supernova kicks have a Maxwellian 
distribution with (Tk > 30kms~^ or so. Our best fit re- 
quires (Tk ~ 15kms~^ though this may be combined to a 
bimodal distribution with a second Maxwellian with (7k2 ~ 
500 kms~^, sufficient to disrupt most systems but to account 
for the high space velocity pulsars. We have also considered 
the distributions of the misalignment between Be star spin 
and orbit that would result from such a kick distribution. It 
is evident that a larger kicks result in larger misalignments. 
We note that the data indicate that the low-velocity kicks 
required to give the current eccentricity distribution might 
not be consistent with the observed misalignments. 

If such a low- velocity kick distribution is ruled out then 
either selection effects must severely limit the observed dis- 
tribution to such an extent that we are only seeing about two 
in seven of the high eccentricity systems or the systems must 
circularise on a timescale similar to their lifetimes. Given 
that this circularisation is not biased towards low periods, 
tides in the stars are not sufficient to be its cause in the 
wide systems. We postulate that a dissipative interaction 
between the neutron star and the decretion disc is a more 
likely mechanism. We have suggested that circularisation of 
Be stars might actually be brought about by an interaction 
between the neutron star companion and the Be star's decre- 
tion disc. This disc is free to grow in size while the neutron 
star is far away. It spends most of its time at apastron. Then 
by periastron the disc may have grown sufficiently that the 
neutron star passes through it and is slowed down. Such an 
interaction can dissipate energy and thus circularise the or- 
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bit. It cannot however change the angular momentum of the 
star nor alter the inclination of the orbit because both these 
processes require transfer of angular momentum as well as 
dissipation of energy. The moment of inertia of the disc is 
smaller than that of the B star and much less than that 
of the orbit. The star-disc interaction can therefore easily 
align the outer parts of the disc with the orbit but any sig- 
nificant change in the angular momentum of the B star or 
of the orbit can only occur on the disc's viscous timescale. 
Meanwhile orbital energy can still be dissipated. Thus the 
misalignment should still contain information on the kick at 
formation even if the eccentricity no longer does. 

Inclinations between the pre- and post-supernova or- 
bits are affected by various factors. Higher kicks lead to 
fewer bound systems but a larger probability of counter- 
alignment. An increase in the mass lost does the same. 
Systems with lower pre-supernova periods are less affected 
and so tend to b e more aligned and more able to survive. 
For the standard iHobbs et al.l (|2005l ) kick distribution Be 
stars ought to be about three times more likely to end up 
with aligned rather than counter aligned orbits but counter- 
aligned orbits should not be uncommon. On the other hand 
if (Tk = 15kms~^ all orbits should be not far from align- 
ment. In order to distinguish between these various possi- 
bilities more information is required about the distribution 
of misalignment angles between the Be star spin axis and 
the orientation of the disc. 
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8 APPENDIX 

Be stars have a decretion disc around them and theory sug- 
gests that Be stars lose their disc to become B stars. If the 
disc reappears the star becomes a Be star again. If we take 
this into consideration there must also be room for the de- 
cretion disc within the Roche lobe of the B star so we insist 
that 



(1 " e)an > C 



Ri 
9il)' 



(52) 



where C is the size of the disc in units of the stellar radius, 
to include a system in ou r analysis. From interferometry 
iGrundstrom fc Gied (|2006l ) find the size of the disc in 7 Cas 
to be Rd =8.1 ± 1.1 R* where R^, i s the r adius of the star 
so C = 8.1 and IGrundstrom et al.l (|2007bt ) find the disc of 
X Per to be about six times larger than the stellar radius so 
C = 6. 

In Fig. [11] we plot the locus defined by equality in equa- 
tion dSg]) against P{ for Ri/g{q) = 6.9 R© for C = 1, 2, 4 
(solid line), 6, 10 and 20 in the case where the initial mass is 
the maximum of AI2 — 25 M© . For each C, the area to the 




1 2 
log,„ (P,/d) 

Figure 11. The solid squares are the Be stars and the triangles 
are the B stars without emission from table [T] The contours show 
e = 1 — ajnin/fln against P{ with C = 1, 2, 4 (solid line), 6, 10 
and 20 from left to right. The initial mass of the exploding star 
is M2 = 25 Mq 



left of the contour is invalid. Only one Be star system re- 
quires C as low as 1.6. This system, 0535-668 is in the Large 
Magellanic Cloud which has low metallicity and hence stars 
have smaller radii. The remainder of Be star systems can be 
accommodated if C = 4. The two B stars which lie to the 
left of the C — 4 line do not have discs because they show no 
emission and so we do not include them in our comparisons. 

One eff'ect of increasing C from to 4 is to reduce the 
number of high eccentricity systems at low period. Our ini- 
tial period distribution biases to low period so this in efi'ect 
reduces P{e) for large e in what follows. Thus if we were not 
to include this our conclusions in Section [5.11 would be even 
stronger. 



